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ABSTRACT 


A computer program, with the code name TRUMP, was developed in 1965 
by the Lawrence Livermore Laboratory to solve a variety of transient 
and steady-state conduction heat flow problems in simple or complicated 
geometric configurations. The purpose of this investigation was to 
adapt "TRUMP" to the Naval Postgraduate School IBM/360 Model 67 com- 
puter system. Several heat conduction problems are solved with the 
adapted version of TRUMP and the results compare closely with the 
analytical solutions to these problems. The example problem inputs 
may be used as guides in preparing input data for future problems. 


Some suggestions are also given for further development of the program. 
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I. INTRODUCTION 


A computer program with the code name TRUMP, was developed in 1965 
by Edwards [Reference H of the Lawrence Livermore Laboratory (LLL) to 
solve a wide variety of non-linear, conduction heat transfer problems. 
TRUMP has been extensively used and continually improved Since its 
initial development and has been successfully adapted to a number of 
computer systems including the CDC/6600 (LLL), IBM/7094 (Bell Labs), 
UNIVAC/1108 (Battelle Northwest, Gulf General Atomics, NBS), and the 
IBM/360 (B.F. Goodrich, NASA Lewis Research Center, Bendix Electro- 
dynamics, Cockerill-Ougee-Providence), to name a few. The objective 
of this thesis project was to adapt TRUMP to the Naval Postgraduate 
School (NPS) IBM/360 Model 67 computer system. 

The B.F. Goodrich Corporation supplied a listing of a version of the 
TRUMP program that they had converted for use in their IBM/360 system. 
Cards were punched from this listing and a version of TRUMP was thus 
secured. 

The B.F. Goodrich version was not complete. In addition to errors 
that resulted from punching 3,350 cards, the program also had an align- 
ment error in one subroutine of the program. Also, several options 
that are available in the complete LLL version of TRUMP were missing 
in the B.F. Goodrich version. 

With a great deal of assistance from the NPS Computer Facility staff, 
the majority of these problems have been corrected and a running version 
of TRUMP is now available at NPS. 

No attempt is made in this thesis to give a complete description of 


the TRUMP program. Rather, limitations of the NPS IBM/360 version as 
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compared to the LLL CDC/6600 version are discussed. In addition, several 
example problems are solved and may be used as guides in preparing input 
data for future problems. For a complete description of the program, 


reference should be made to the TRUMP users manual [Reference iJ 5 
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If. DESCRIPTION OF THE PROGRAM 


A. GENERAL 

The program solves the general, non-linear, parabolic partial dif- 
ferential heat equation. Transient and steady-state heat conduction 
problems in multidimensions may be solved. Initial conditions, material 
properties, source and sink strength and boundary conditions may vary 
with time, spatial position and/or temperature. 

In order to solve problems with TRUMP, the user specifies input 
data in "BLOCK" form to 12 input data blocks. A description of each 
of these blocks and their use will be given in Section B 2. 

Following the definition of the problem a model must be constructed. 
Complex geometric shapes are divided into regions and if there is sym- 
metry, the symmetry planes are represented by insulated surfaces in 
order to simplify the model. Regions may have different materials with 
specified properties. 

The regions may be subdivided into volume elements or nodes having 
arbitrarily any size and shape. Nodes or volume elements must have a 
representative nodal point whose location may be anywhere in the node 
or on the surface of the node. However, in transient problems, for 
maximum accuracy, node shapes and nodal point locations are chosen so 
that the lines joining the nodal points of connected nodes are perpen- 
dicularly bisected by the connected area. In steady-state problems, 
due to the lack of heat capacity effects, the nodal points may be lo- 
cated anywhere in the node without loss of accuracy. A variety of 
boundary conditions may be given to the surface nodes. Constant temper- 


ature, variable temperature, constant heat flux, variable heat flux, 
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forced convection conditions, free convection conditions, radiation 
boundary conditions or combinations of the above may be treated. 

Instead of having very fine zoning of the regions, rather crude 
zoning will give the solution with less effort. Later, the solution 
may be obtained for a finer subdivision of the system. Also, with 
the crude zoning one can easily understand which regions need finer 
zoning. 

The accuracy of the given problem depends on how well the problem 
can be modeled, and the fineness of the spatial subdivisions and time 
steps used in the calculations. 

Allowable problem size limits for the NPS IBM/360 version of TRUMP 


are summarized in Table l. 
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TABLE I 


Problem Size Limits 


Item Data Block Parameter Size 
Materials 2 _ My HS 
Reactants 8) M3 5) 
Nodes 4 My 355 
Internal thermal connections 5 Ms 950 
External thermal connections 6 Me 60 
Boundary nodes 7 M7 20 
Heat generation tables 8 Mg 5 
Initial conditions 9 My, 35) 
Mass flow connections 10 Mio 50 
Temp. versus time plot nodes ite M1 10 
Temp. versus plot times > M) 1 
Remotely dependent properties WW Mj 9 75 
Table lengths DS, WOR ORO Mg 152: 
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In general, any consistent set of units may be used for the input 
data; output quantities are then given in the same unit system. 

Each material, system node, and boundary node can be identified 
by a non-zero integer number fitting within 5-column fields on the 
data cards. These identifications can be in any sequence. [In order 
to reduce the number of input cards arithmetic sequencing is suggested 
so that the nodes having identical descriptions can be submitted as 
input on only one card. 

For non-symmetric, symmetric with respect to an axis, or to a 
center, the nodal volumes and the thermal connection areas can be cal- 
culated by a control value in the input data. 

Initial conditions may be specified in BLOCK 1 of the input data. 
These values will be assigned to all nodes described in BLOCK 4, unless 
other values are specified for individual nodes in BLOCK 9. 

Each internal thermal connection between nodes must be described 
in the input data by specifying the two node identification numbers, 
two connector lengths, and two interface dimensional factors. The 
selection of connector lengths depends on the location of the nodal 
point within each node. These locations may be arbitrary, but in 
transient problems should usually be at the geometric centers of the 
nodes. If the amount of the heat flow depends on surface temperatures, 
then they can be located on the surface. 

A system that exchanges heat with its surrounding must have external 
temperatures specified along its boundary. In the input data, boundary 
node identification numbers and their temperatures, which may be con- 
stant or a function of time, must be given. The term "boundary node" is 


identical in meaning to an external temperature. 
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Surface temperatures can be specified by using boundary nodes with 
specified temperatures and connecting the surface nodes to a boundary 
node with a very large convection coefficient. In order to reach a 
high degree of accuracy, surface nodes should be selected as eitler 
zero-volume nodes, or very thin nodes. When a large number of surface 
nodes are connected to the same boundary node and the heat transfer 
coefficients are not tabulated, the number of external connections can 
be reduced by replacing them with connections to a single zero-volume 
node that is connected to the boundary node with a very large convection 
eoefficient. 

In order to generate steady-state solutions, TRUMP must first go 
through a set of transient calculations. The steady-state solution is 
obtained by allowing enough problem time to take the transient to its 
conclusion. However, if only steady-state results are desired, the 
control on the average temperature change during a time step (TVARY) 
can be made very large. This allows the program to neglect transient 
stability limitations on its choice of time steps resulting in a mini- 


mum expenditure of computer time on transient calculations. 


B. USE OF COMPUTER 
1. General 
The data deck must include a problem name card, any number of 
block number cards with their input cards, and a data end card. The 
first card of each datadeck must be a problem name card. The symbol "*" 


must be in column 1, any identification of the problem in columns 2 


through 71 and 73 through 80. 
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In the program there are 12 allowed input data blocks, with 
block numbers from 1 through 12. Each block must begin with a block 
number card, and except for BLOCKS 1 and 11, must end with a blank 
Bea. The block number card must begin with the word "BLOCK" in 
columns 1 through 5, and the block number, from O01 through 12, in 
columns 6 and 7. 
Data BLOCK 1 must follow the problem name card; all other data 
blocks may be placed in any order. 
The last card of each data deck must be a data end card, with 
fhe word "ENDED". in columns 1 through 5, and "-1L" in columns 6 and 7. 
Additional cards containing comments without '*" in column 1 
may precede the name card and will appear on the printout. Any desired 
block description in columns 9 through 80 may be made on the block 
number card, This will also appear on the printout. 
2 block, Lien Descriptions 
Basically, one should refer to the TRUMP users manual [ Reference 1 
for format details and the original descriptions of the input values. 
In this section the main difference between the IBM/360, Model 67 version 
of TRUMP and the LLL, CDC/6600 version are discussed. 
in, BLOCK I: 
MCYC Maximum allowed number of time steps. If zero, neg- 
ative, or unspecified will be 30,000. 
MSEC Maximum allowed machine time, in seconds. If zero, 
negative, or unspecified will be 30,000. 
NPUNCH Indicates that the final values of temperatures of nodes, 
weight fractions of reactants, and constant heat gener- 


ation rates of nodes will be punched. Since punched 


ES) 


athe ,edeesd atab teqnd ieealls B1 asc erm eeigeny e484 ut ; 

dcoafd @ €2)0 algec tame M6019 asa2 tr sawenl? 2D movil ‘epedsee avebd 

son f4 a dihv hue Gewese ,4f bos 2 GOON te) sqnone bem .e2as ¥ 

of *SD00E" “beow edls dedw'nigad teu brad sedens Sadid att eg 

nl [Sf diguesds 10) wort ,rodawsr dookd oda bre ,¢ figueads J 

* tee 6 

wieb sari80 ile Lbr8> amen milder sf3 wolfe] scum f WO1d ast 

4.2. ,tebto qm ae) bemiq ed Yan 

stv , 6300 -ba9 nikb » od te0m dosh 230 isis 0 bees geal ofF 

© Gun 6 anmplos g4°%F-" bon @ dauonts 1 anwios nt "Osc t 

~ f P emutoo axe Sundaae etnemnon. grtninsion sins Lasohaabted 
7 teerdeeh yuh ee Tae {fiw bre bse omen ais 4 
“daotd sdd°nis: Oban 2 enn: 08 iiquércits ©: aime lon nk wehaqhebe 


- 


ae. 

ae © “.dneaiag of3 no seegqs. cele LLiw ent .t 
- z id 7 ; all | ORE ee - nosIgrr{eg = .J 9D, 

: | a ee inotand 


y utes sagih-oAs 20: smbbigtiro set fonikgizo add ban wELdeb' 126 
Madexew (4 ladah DAC WEE pets aewsed, conessth 2b: chem odd ) 
-boseuneds isa ralexev OO80\C2 ~LLT eis wa. 
ee. | HOI & 7 
“yao .oyss af. .eqe34 ‘east to 1sdmm bewolls auaixes. 
Se od itiv beltiveqenu 2 avidin 
eek niin a ants gabdsse havelle aumixad ORE 
. 000, of od thaw bold} seqanv 70 javiseaee 
ie lo @etatassgm6s 26 ebufer aunt ott tad” ansaolbal OM 
ahs om ted -tnasen6> bow ,sanedoney le encison7? Jilgiav 


bvlseoy aoatt | badshug'ed Dilv eaberae go3n2 eolte - 


21. 


a =< 7) 


output is not obtained in this version of TRUMP, the 
value of NPUNCH should not be used. 

IRITE Number of time steps between data output on unit ITAPE 
of problem time and table of node numbers and node 
temperatures, in addition to output of first, second, 
and last time steps. Only used in TRUMP/6600/DS version. 
This) control is not in use in this version. 

ITAPE Output unit on which output data will be written as 
described in IRITE, This input value is not in use 
for IBM/360 Model 67 version of TRUMP. 

TIMEP Problem time interval between data output, in addition 
to output on first, second, and last time steps. This 
input control value is also not in use for this version 
of TRUMP. 

In addition to the above differences, the data transfer option 

Bsonot in use in this versig of TRUMP. 
3. Numerical Output 
The amount of output data is controlled by the value of KDATA, 
However, for every problem, all input data and several other quantities 
such as node volumes and connection areas are written out. The results 
of the first, second, and last time steps are always written out. In 
addition, output data are written out at problem time step intervals 
of IPRINT. For the notations used on the printouts one should refer 
to TRUMP users manual [Reference 1 : 
4, Evaluation of Results 
Any desired end to the problem may be made with the specified 


input values in BLOCK 1. These may be maximum problem time, TIMAX, 
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maximum temperature, TMAX, minimum temperature, TMIN, maximum number 
of time steps, MCYC, or maximum number of seconds of machine time, 

MSEC. On the printouts, the value of the BOR rence sentinel, KWIT, 
will appear. Table II shows the values of KWIT, and the descriptions 


of these values. 
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TABLE II 


Values of KWIT 


KWITt Cause of Problem Ending 
i The problem time, SUMTIM, reached TIMAX 
Z A temperature exceeding TMAX 
3 A temperature less than TMIN 
4 Steady-state 
5D Data 
6 Relates to the tape usage 
7 The number of time steps reached MCYC 
8 Machine time reached MSEC 
9 BLOCK 2 or BLOCK 4 missing 
10 Convergence failure 
11 Items exceeded the size limits 
2 Table lengths exceeded the size 
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DESCRIPTIONS OF SUBROUTINES 


THERM 


TALLY 


CHEM 


SPECK 


FINK 


GEN 


SURE 


FLOW 


SEEK1 


SEEK2 


Subroutine for material properties and node descrip- 
tions. **BLOCKS 23°45 Vand 12 are included: 

Subroutine for initialization, totaling, and checking 
results, converting regular to special nodes, testing 
for ending the problem. BLOCKS 1 and 9 are included. 
Subroutine for chemical reaction. BLOCK 3 is included. 
Subroutine for heat flow between special nodes and 
other nodes or boundary nodes. [Iterates to solve 

set of implicit difference equations when special 
nodes connected to each other. 

Subroutine for heat flow between nodes by conduction, 
convection, and radiation. Treats all nodes as regu- 
lar nodes. BLOCK 5 is included. 

Subroutine for internal heat generation. BLOCK 8 

is included. 

Subroutine for heat flow between surface nodes, bound- 
ary nodes by radiation, free and forced convection 
treats all surface nodes as regular nodes. BLOCKS 6 
and 7 are included. 

Subroutine for mass flow between nodes. Treats all 
nodes as regular nodes. BLOCK 10 is included. 
Subroutine for finding an identification number of 

a material, reactant, node or boundary node in a 
block item list. 

Subroutine for finding the identification numbers 


of a pair of nodes ina: block item dist. 
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MAINT 


TLE 


CLOCK, CLOCKL 


PLOT 


PATCH 


REFER 


SRCON 


Subroutines for initialization before first time 
step. 

Subroutine for making CRT plots and tables of temp- 
erature versus node location, temperature versus 
time, and node location versus time. BLOCK 11 is 
included. 

Subroutine for converting a number with a decimal 
point to floating number, otherwise substituting 

a specified floating point number. 

Subroutine for finding array subscripts of materials, 
reactants, nodes or boundary nodes referred to by 
identification number in another data block, and 
writing out diagnostic statements whenever one can 
not be found. 

Subroutine which is related to PATCH. Converts a 
number having decimal point to a floating point 


number. 
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IIIT. EXAMPLE PROBLEMS 


A. SLAB-MELT PROBLEM 

A semi-infinite slab (-co€ x Coo, - cody Coco, 0 gz <coo), 
has thermal conductivity 1.0 cal/sec-°C, density 1.0 g/com>, heat capacity 
LO cal/sec-°C, a melting point of 50.0°C, and a latent heat of fusion 
of 25.0 cal/g. The initial temperature is 100.0°C, and for times E> OF 
the surface, z = 0, is maintained at 0.0°C. TRUMP is used to calculate 
the temperature distribution in the slab, and the motion, of the solid- 
liquid interface, Zint, for times up to 200 sec. The results are com- 
pared with the exact solution. 


1. Analytic Solution 


From Reference 3 analytic results are: 


iL 
fpy eel 2eaient (peg) for) 0 GE 50eG 
ale 
Z 
te 0 4 0. Wiers (=) | for) 50 Ge 100°C 


att 
Zint = 0.835:{t for T = 50°C 


2. Trump Solution 

The input data for this problem are shown in Figure 1. In order 
to have output at every 10 cycles, in addition to output on the first, 
second, and last time steps, IPRINT is selected as 10 in BLOCK 1. For 
temperature accuracy, TVARY is 1.0°C and the problem is ended at 200 sec 
by setting TIMAX equal to 200.0. TONE is 100.0 since the initial tem- 
perature of the slab is 100.0°C. Homey since the zero-volume sur- 
face node, node -1, is at 0.0°C for t > 0 its temperature is specified 


in BLOCK 9 by setting the value of TT equal to 0.0. In BLOCK 2, the 
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material name, AMAT, and material identification number are arbitrarily 
selected at SMELT and 50, respectively. Other slab thermal properties 
such as density, heat capacity, thermal conductivity, temperature at 
which the latent heat is released and the latent heat of fusion are also 
specified in BLOCK 2. In order to give node descriptions in BLOCK 4 


the slab is divided into nodes as shown below. 


Gila 


A total field depth of 100.0 cm is chosen so that at 200 sec no 
appreciable temperature change will take place at the maximum value 

of Z. In this figure, the node identification numbers, and the node 
thicknesses are shown. The nodal points of all nodes are at the node 
centers. In BLOCK 5, internal connection node numbers, lengths of 

heat conduction paths from the nodal points to the connected interfaces, 
and the area of the connected interfaces are defined. In BLOCK 6, 

a zero-volume surface node, node -1, is connected to a boundary node, 


node 2001, with a heat transfer coefficient of 1.0 x 108cal/sec-cm2-°C, 
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and the external surface area of the surface node is selected as unity. 
In BLOCK 7, a constant boundary node temperature is defined as 0.0°C. 


Computer running time for this problem was 65.52 sec. 
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FIG, 1 INPUT DATA FOR SLAB MELT 
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‘TEMPERATURE DISTRIBUTION IN THE SLAB-MELT PROBLEM 
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TABLE III 


SOLID-LIQUID INTERFACE IN THE SLAB MELT 
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B. FINITE-SLAB PROBLEM (suggested by Professor Paul Pucci) 
A Large, lO vineh“thiek "steel plate, initially at "80.0°F, is 

suddenly submerged in a fluid at temperature 680.0°F. TRUMP is used 

to calculate the transient temperature distribution in the plate assuming 
a constant surface heat transfer coefficient of 1248.0 BTU/hr-f£t2-°F. 
The thermophysical properties of the steel were taken as: thermal con- 
ductivity 26.0 BTU/hr-ft-°F, density 487.0 1bm/ft?, and specific heat 
0.133 BIU/1bm-°F. 


fT. Analytic Solution 


i a oO Sind,L - (nL) 2 Fo 
e Cos (AnL=) 


mee 
ite POS ApL + Sind,LCos AL 


hL 
Applying the Biot number, Bi = = 2.0) to the e1genivaliue 


AnL 
problem of Cot (AnL) = ; the two eigen values AL = 1.0769 and 
Bi 
ML = 3.6436 can be found. After substituting the appropriate value of 


—_—— 


the Fourier number, Fo = io Ss 230.4t in the above series solution and 
considering only the first two terms, the result becomes: 


jt) + 680 - 1200|0.58923153 qr Ous2e Cos (1.07698) - 6.11es6085e cos (3.6436) 


Zea Lrump Solution 

The input data for this problem are shown in Figure 3. For unit 
consistency, lengths in inches are converted to feet and applied to the 
problem. Since the plate has symmetry with respect to the mid-plane, 
this plane is taken as an insulated plane in the TRUMP solution. So, 
in BLOCK 1 symmetry type indicator, KD, is applied as 1 for the non- 
symmetric case, and the total thickness of the steel plate is considered 
as 0.5 inches. Temperature scale indicator, KT, is selected as 3 


because of the Fahrenheit scale. For temperature accuracy, the desired 
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maximum temperature change in each time step, TVARY, is selected as 

1.0°F. Maximum problem time, TIMAX, is given as 1.0 hr. TONE is set 
equal to 80.0 making the initial temperature of all nodes equal to 80.0°F. 
In BLOCK 2, material name and material number are given as, FSLAB and 

50, respectively. Also, the steel properties such as density, specific 
heat, thermal conductivity are defined. In BLOCK 4 node numbers, mate- 
rials of nodes and their dimensions are given with respect to the equally 
spaced 10 nodes for the 0.5 inches of thickness of the slab. In BLOCK 6, 
the zero-volume surface node, node -1, is connected to a boundary node, 
node 2001, with a surface heat transfer coefficient of 1248.0 BTU/hr-ft2-F. 
In BLOCK 7, a boundary node temperature of 680.0°F is given. The com- 
puter running time for this problem was 24.46 sec. 

In order to see the effect of finer zoning on the results, the 
same problem was also solved using 20 and 40 nodes. However, the results 
were all within 1°F accuracy and not much improvement was obtained with 
the finer grid sizes. Computer running time, on the other hand, was 
@mereased from 24.46 sec. to 53.76 sec. and to 172.64 sec. for*l10, 20, 


and 40 nodes, respectively. 
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FIG. 4 TEMPERATURE DISTRIBUTION IN THE FINITE SLAB PROBLEM 
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C. SEMI-INFINITE SLAB PROBLEM 

A semi-infinite slab is initially at a uniform temperature of 100.0°C. 
At time t = 0, the surface temperature is suddenly lowered and main- 
tained at a temperature of 0.0°C. The slab has a thermal conductivity 
o£ 1-0\ cal/sec cm=-°C, density of 1.0 g/cm, and a heat capacity of 
1.0 cal/sec-°C. The temperature distibution in the slab, and the total 
heat flow from the slab are calculated as a function of time for times 
up to 200 sec using TRUMP, and the results are compared with the ana- 


iyewcal solutions. 


1. Analytical Solution 
x2 


| 


rox, ct) = 200 ert 
At 


= 
Qe = -112.84 | t 


2. Trump Solution 


The input data form is shown in Figure 6. In BLOCK 1, problem 
controls such as problem time, TIMAX = 200.00 sec, maximum temperature 
change in each time step, TVARY = 1.0. G. and the initial temperature of 
alll modes, TONE = 100.0°C are given. In BLOCKS 2, 4,5, 6, 7 and 9, the 
thermal properties, zoning procedure, identification numbers, internal 
and external connection descriptions, boundary node temperatures, and 
initial conditions are defined identically as in the slab-melt problem. 


Computer running time for this problem was 27.05 sec. 
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D. CYLINDER PROBLEM 

An infinitely long cylinder of unit radius, with a unit surface 
temperature, is embedded in a medium with a zero initial temperature. 
TRUMP is used to calculate the temperature distribution in the medium. 
All thermal properties for the medium have unit values. Results are 
compared with an analytical solution. 

L.. rump. Solution 

The input data are shown in Figure 8. Because of the geomet- 

rical shape of cylinder, the symmetry type indicator, KD, in BLOCK l 
is selected as 2. Maximum problem time, TIMAX, is made 50, and the 
medium initial temperature, TONE, is set to 0.0. Also, for temperature 
aceuracy, TVARY, is selected as 0.05. In BEOCK 2, unit properties of 
the medium are defined. Nodal point locations are chosen so as to be 
able to compare the TRUMP solution with the given table in Reference 5 


rather easily. The zoning of the medium is shown in the figure below. 
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In BLOCK 4, the arithmetic mean average radii, DRAD, are used to obtain 
correct nodal volumes for cylindrical shaped nodes. Actual node inter- 
faces are used as interface radii, DRAD, in BLOCK 5, to obtain the areas 
of the connected interfaces. In BLOCK 6, zero-volume surface node, 

node -1, is connected to the boundary node with a heat transfer coeffi- 


elent of 1.0 x 108, and geometric factors for external surface node 


area are defined. Computer running time for this problem was 15.46 sec. 
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FIG, 8 INPUT DATA FOR CYLINDER 
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TABLE V 


TEMPERATURE DISTRIBUTION FOR CYLINDER (TIME = 50) 


R T (Trump) T (Analytic) 
0 1BELO) 10 
io 0.96% 0.963 
LS 0.3897 0.898 
ee) 0.842 0.843 
Ay’ 0.794 0.794 
low) Ocioe Os/5L 
310 0.574 0.574 
D0 0.382 0.381 
7.0 OsZ6r 0.260 
a0 O79 Oy, 
10.0 0.146 0.146 
20.0 0.016 0.016 
30.0 0.001 0.001 
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E. SPHERE PROBLEM 

A sphere with a unit radius and a unit surface temperature is 
embedded in a medium with a zero initial temperature. The temperature 
distribution in the medium is calculated using TRUMP assuming all the 
thermal properties have unit values. Results are compared with the ana- 


ivtereal solution. 


1. Analytic Solution 
t 
iG) = enfeG— =)) 


2. Trump Solution 


The input data from the cylinder problem can be used to solve 
this problem if the symmetry type indicator, KD, in BLOCK 1 is changed 
from 2 to 3 to correspond to the centrisymmetric geometry. All other 
input conditions are identical with the cylinder problem. The computer 


time used for this problem was 11.61 sec. 
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TABLE VI 


TEMPERATURE DISTRIBUTION FOR SPHERE (TIME = 50) 


R T_ (Trump) T (Analytic) 
0 1.0 1.00 
ieee 0,89 0.90 
MGS) 0.75 0.75 
ee, 0.64 0.64 
pe 7 0.56 0256 
9 0.49 0.49 
2.0 0.28 0.28 
Bie) Ore 0.14 
iO 0.08 0.08 
9.0 0.05 0.05 
10;.'0 0.04 0.04 
20.0 0.00 0.00 
30.0 0.00 0.00 
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F. TWO DIMENSIONAL PLATE PROBLEM 

A rectangular plate, 40 x 20 cm, has a thermal conductivity of 
400 cal/cm-sec-°C, a density of 1.0 elem - and a heat capacity of 
1.0 cal/e-°C oe is) initially at 0).0°C.- “Three sides.of the plate are 
maintained at 0.0°C. The upper side of the plate has a temperature 
distribution that varies sinusoidally with x, starting at 0°C, reaching 
a maximum of 100°C at the ante and decreasing to O°C at the outer 
edge (see figure below). The steady-state temperature distribution is 
calculated in the plate using TRUMP. Results are compared with the 


analytical solution. 


Tatoo sin AG) 


- 
i} 
O7 
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| | T=0°C 
40cm 


ils Analytic Solution 


f TX 
Ge, = 100 Sin (— 
(x,y) nig 
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2. Trump Solution 


The input data are shown in Figure 9. Because the boundary 
conditions of the plate and the plate itself are symmetric with respect 
to the plane at x = 20 cm, this symmetry plane is assumed to be an in- 
sulated surface. The zoning of the plate is shown in Figure 10. [In 
BLOCK 1 of the input data, IEPRINT, is selected as 9999, so that print- 
outs are obtained only for the first, second, and last time steps. To 
obtain a rapid conclusion to the transient problem, TVARY, is selected 
as 100.0°C. In BLOCK 4, instead of using a very large number of surface 
nodes for sides x = 40 cm and y = 0, surface node -1 and surface node -4 
are used, respectively, in order to reduce the number of the input cards, 
and in BLOCK 6 they are connected to the boundary nodes 2001 and 2004. 
In BLOCK 7, the boundary node temperatures, and in BLOCK 9, the initial 
temperatures of the surface nodes are given. The computer running time 


for this problem was 28.71 sec. 
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IV. CONCLUSIONS 


The adaptation of the computer program "TRUMP" to the Naval Post- 
graduate School IBM/360 Model 67 has been achieved, with success. 
Analytical temperature distributions of various transient and steady- 
state heat conduction problems compared very well with the results 
obtained from TRUMP. The complete adaptation has not yet been estab- 
lished, however. At present, punched output and plotted graphs have 
not been obtained. The computer solutions of the example problems did 
not require the use of subroutines CHEM, GEN, and FLOW. Thus, the 
adaptation of these subroutines is yet to be evaluated. Since, the 
accuracy of the solutions are dependent on the zoning, this effort can 
be reduced if an existing program "FED: A computer program to generate 
geometric input for the heat-transfer code TRUMP" by Dale A. Schauer 
‘Reference 4 is adapted to the IBM/360 Model 67 version of TRUMP. 

The example problems presented herein will provide a useful reference 


for preparation of input data for future problems. 
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